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Abstract 

The vibrational properties (phonons) of a one-dimensional periodic phase 

of atoms physisorbed in the external groove of the carbon nanotube bundle 

are studied. Analytical expressions for the phonon dispersion relations are 

derived. The derived expressions are applied to Xe, Kr and Ar adsorbates. 

The specific heat pertaining to dense phases of these adsorbates is calculated. 
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I. INTRODUCTION 



Adsorption of gases in nanotube materials has been investigated extensively in recent 
years, both from the experimental^"^ and theoretical^"^^ standpoints. A very special ge- 
ometry of the nanotube materials provides a possibility for realization of unique adsorbate 
phases with reduced effective dimensionality^'^^'^^'^^. The behavior of adsorbates depends 
on whether they are adsorbed within the tubes^, in the interstitial channels^, or in the 
grooves. For samples consisting of nanotubes with closed ends, the adsorption in the 
tube interior is not possible, and for sufficiently large adsorbates, the adsorption potential 
in the interstitial channel is purely repulsive.^ Thus, for the case of large atoms adsorbed 
in the material made of closed-end carbon nanotubes, the only positions available for gas 
adsorption are the grooves of the bundles of carbon nanotubes. 

In this paper we consider a dense, periodic phase of noble gas atoms physisorbed in the 
groove positions of a nanotube bundle. This phase can be visualized as a one-dimensional 
chain of atoms arranged in a periodic array within a groove (see Fig. 1). While the single 
particle properties of quantum adsorbates in grooves have been considered in Ref. 10, in 
this article we concentrate on the many-particle excitations (phonons) of the dense, periodic 
adsorbate phase, that crucially depend on the interactions between adsorbates. 

Recently, measurements of the specific heat of atoms adsorbed in nanotube materials 
have been reported^'^. In this respect, it is important to theoretically consider the specific 
heat of a dense phase and its signature in the overall specific heat of the sample (nanotube 
material -|- adsorbates). 

The outline of the article is as follows. In Sec. II, we present a simple lattice dynamics 
approach to the vibrations of adsorbates in a groove. We show how the phonon frequencies 
are related to the interadsorbate interaction potential and the interaction of an adsorbate 
with the surrounding nanotube medium. In Sec. Ill, we apply the model from Sec. II to 
adsorption of Xe, Kr and Ar atoms and we calculate the characteristic vibrational frequencies 
pertaining to these adsorbates. In Sec. IV, we use this information to calculate the specific 
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heat of dense adsorbate phases. Sec. V summarizes the main results. 

II. VIBRATIONS OF DENSE PHASES OF ADSORBATES IN THE GROOVES 

To study the vibrations of the adsorbates in a groove we introduce a simphfying assump- 
tion of rigid substrate i.e. we do not consider phonons of the nanotube material. In this 
case, the total potential energy $ of the adsorbate system is given as 



Here, v{ri — r;/) represents the interaction between the adsorbates at r; and r;/ positions. 
The indices I and I' denote the adsorbed atoms, and V{ri) represents the interaction of the 
l-th adsorbate with the surrounding nanotube medium (substrate). We assume that in the 
absence of vibrations, the adsorbates are positioned in an infinitely long one-dimensional 
lattice (with the lattice parameter a, see Fig. 1) within a groove. We also neglect all the 
interactions between the adsorbates which are not in the same groove. This is an excellent 
approximation, due to a large diameter of carbon nanotubes. 

If the interaction with the substrate were vanishingly small, the problem would reduce to 
three dimensional vibrations of an unsupported chain of atoms. This problem can be easily 
solved analytically as shown e.g. in Ref. 15. Even with the external field V, the problem can 
be analytically solved, following the lattice dynamics approach. The eigen-frequencies 
of the adsorbate system can be found from the matrix equation 



where M is the adsorbate mass, a;(g) is the vibrational frequency, q is the phonon wave 
vector parallel to the groove and e{q) is the 3D phonon polarization vector. Note that the 
phonon wave vector g is a one-dimensional quantity since the adsorbate lattice exhibits a 
ID translational invariance. The polarization vector e(g) is a three-dimensional quantity 
since the displacements of adsorbate atoms in all three spatial directions are allowed. If we 
take into account only pair interactions v between the nearest neighboring adsorbates and 




(1) 



[V{q)-Mu;{qf]e{q)^0, 



(2) 
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neglect the discrete nature of the tubes surrounding a groove by effectively "smearing" the 
carbon atoms along the tube surface^'^, the 3x3 dynamical matrix, V{q) is given as 



2/3[l-cos(ga)]+/, 



XX 



xy 



xz 



xy 



2a[l — cos(ga)] + /j 



yy 



(3) 



xz 



yz 



2a[l - cos(ga)] + /, 



zz 



The X-axis is oriented along the chain of adsorbates, y-Sixis passes through the centers of the 
two tubes surrounding a groove, and the z-axis in perpendicular to both x and y axes and 
points outward from the bundle, as denoted in Fig. 1. The external force constants /^^,^ are 
given by 



where the subscript zero signifies that the derivatives should be taken at the equilibrium 
position of the adsorbate. As the tubes surrounding a groove are assumed to be smooth, 
the external potential, V{ri) does not depend on the adsorbate x-coordinate and there is 
no external (substrate induced) restoring force associated with the adsorbate displacements 
in x-direction. That is why all /^^i^ constants with = x or u = x vanish in Eq. (3). 
Additionally, due to the symmetry of the external potential with respect to the displacements 
parallel to the line joining the centers of the two tubes surrounding a groove (|/-direction, 
see Fig. 1), the constant fyz also vanishes. 

The force constants associated with the adsorbate-adsorbate interactions [a and P in Eq. 
(3)] are given as 



where r is the coordinate of relative distance between the neighboring adsorbates. The 
phonon dispersion relations, uj{q) for —n/a < q < Ti/a, obtained as the solutions to Eq. (2) 
with /^^ = f^y = ^ fy^^O are given as 
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(4) 




(5) 



UJT2{q) 




(6) 



where the subscripts L,T1,T2 denote the longitudinal (L) and the two transverse phonon 
branches (T1,T2). The L, Tl and T2 modes are polarized exclusively in x,y, and z di- 
rections, respectively. Similar equations have been derived in Ref. 12. The authors of this 
reference considered an isotropic oscillator model, V{x, y, z) — f[{y — yQy-\-{z — ZQY\/2. The 
two transverse modes are degenerate in this approximation, which however is not always the 
case, because generally fyy ^ f^z, as will be shown in the next section. 

A brief comment concerning the neglect of corrugation of the holding potential, i.e. 
discreteness of the nanotube, is in order here. If the adsorbate phase is commensurate with 
the holding potential provided by the substrate, the longitudinal mode will exhibit a zone 
center gap which can be related to the magnitude of the corrugation. This effect has been 
experimentally confirmed for Xe overlayers on Cu(lll) surface^^'^^ and should in principle 
be observable for all commensurate systems. 



The ingredients necessary for the calculation of phonon modes are obviously the force 
constants a,/3 and fyy, fzz- The first two can be calculated from the presumably known 
interadsorbate interaction potential, while the /-force constants can be calculated from the 
asdsorbate-substrate interaction potential. 

The adsorbate-substrate interaction can be calculated as a superposition of the 
adsorbate-nanotube interactions. It is sufficient to consider only the two nanotubes sur- 
rounding the groove, since the interaction of adsorbate atom with other nanotubes is van- 
ishingly small due to a large diameter of the carbon nanotubes typically encountered in 
experiments.^^ Neglecting the discrete nature of the carbon nanotube, the adsorbate-single 
wall nanotube interaction can be written as'^'^ 



III. VIBRATIONS OF Xe, Kr AND Ar ADSORBATES 
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V{p) = 3n9ea^ 



"21 




32 
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(7) 



where e and a are the energy and range parameter of the effective adsorbate-carbon site 
interaction which is assumed to be of a Lennard- Jones form. The variable p denotes the 
distance of the adsorbate atom from the axis of the tube, 9 is the effective coverage of C 
atoms on the tube surface (^=0.38 1/A^), R is the tube radius and the variable rj is defined 
as T] — R/ p when p > R. The function Mn{rj) is defined as 

d(f) 

"^-^"^ ^ Jo (l + ,^-2,cos0)"/^- 

The external potential can be constructed as a sum of the adsorbate interactions with two 
tubes surrounding a groove. Prom the thus constructed potential, the relevant external force 
constants can be obtained using Eq. (4). 

To illustrate the shape of the holding external potential, we plot in Fig. 2 the potential for 
Xe adsorbate along the y and z directions. In this figure and in all subsequent calculations, 
we consider single wall carbon nanotubes with the diameter of 13.8 A , arranged within a 
bundle in a triangular lattice with the lattice constant of 17 A. The origin of the coordinate 
system is chosen in such a way that the centers of the two tubes surrounding a groove are 
positioned at (7/ = A, 2; = A) and (y = 17.0 A, 2; = A). The parameters of the 
Xe-C site potential are obtained from the so-called combination rules, as explained in Ref. 
7. We adopt all the relevant parameters of adsorbate-adsorbate and adsorbate-substrate 
interactions from Table I of Ref. 7. Since the two curves displayed in panels (a) and (b) of 
Fig. 2 pass through the absolute minimum of the external potential {y = 8.5 A, z = 6.38 A), 
the external force constants can be obtained as second derivatives of the potential curves 
presented in Fig. 2 at their minimum positions. Note that the potential is significantly more 
anharmonic in the z than in the y direction which means that the harmonic approximation 
to the potential is less satisfactory in the 2;-direction as can be seen from inspection of thin 
dashed lines in Fig. 2 which represent the harmonic approximation. It should be mentioned 
that the minimum value of the potential we obtain is consistent with the one reported in Ref. 
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7. However, the binding energy of a single Xe adsorbate in this potential is not consistent 
with the experimental result^ (-282 ±11 meV) which is lower than the absolute minimum of 
the potential we obtain here (-225 meV). A possible small difference in the tube separations 
at the surface of the bundle with respect to the corresponding separation in the bundle 
interior may be responsible for this discrepancy. In addition, as shown and discussed in Ref. 
10, small changes in the adsorbate-carbon site effective potential may cause relatively large 
change in the adsorbate binding energy. 

The total adsorbate- adsorbate interaction encompasses direct adsorbate-adsorbate in- 
teraction and the indirect interaction mediated by the polarizable substrate material and 
other adsorbates.^^'^^ This interaction is generally different from the corresponding interac- 
tion between the atoms in the gas phase. The most critical adsorbates to be considered in 
the remainder of this paper are Xe adatoms which have the highest polarizability of all the 
adsorbates we shall consider. While the substrate mediated forces can be calculated for a 
planar substrate material, it is more difficult to calculate them and assess their importance 
for the substrate of present interest, i.e. a bundle of carbon nanotubes. This issue has been 
discussed in the case of interstitial adsorption in Ref. 20. 

Even though the many body forces should play a role in dynamics of Xe monolayers 
on crystal surfaces^^, it is often found that the gas phase Xe-Xe potentials provide very 
reliable phonon dispersions^^'^®. This is in part due to the fact that phonon dispersions are 
not determined by the potential itself, but by its derivatives which may be less influenced 
by many-body effects. The only exception to this "rule" seems to be the somewhat puz- 
zling dynamics of Xe over layers on Cu surf aces. ^^"^^ In all the calculations to be presented, 
we shall assume that the adsorbate-adsorbate interactions can be well represented by the 
corresponding gas phase potentials. 

We shall model the adsorbate-adsorbate interaction also by the Lennard- Jones potential, 
4^33 [('^gg/'")''^^ ~ ('^ss/'")^]) where r is the interadsorbate separation, and egg and agg are the 
energy and range parameters of the potential, respectively. The equilibrium positions of the 
adsorbates can be found from the static equilibrium condition, i.e. by minimizing expression 
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(1) with respect to the lattice parameter a. Following this procedure, we find that the ID 
lattice parameter is given by 

a = a„„7r ( I ^ 1.1193(7o„. (9) 

V 6756757 ^ ^ 

The force constants a and P [Eq. (5)] are then given as 

a = -0.1665%- 

(5 = 60.615%. (10) 

Force constant a determines the dispersion of the transverse modes. It is negative due 
to the fact that the lattice parameter a is smaller than the equilibrium separation of the 
binary Lennard- Jones potential (1.122 Ggg). Note also that \a\ <^ \/3\, which means that the 
transverse modes have negligible dispersion in comparison with the longitudinal mode. 

If the adsorbates form a chain that is commensurate with the substrate corrugation, the 
periodicity of the one-dimensional chain of atoms may be different from the one obtained 
here by neglecting the substrate corrugation. In this case, all the force constants may change. 
On the other hand, for the incommensurate chain, our model should provide reliable results. 

In Table I, we summarize the lattice parameters and the force constants needed for 
calculation of the phonons in the specific cases of Xe, Kr and Ar atoms adsorbed in a 
groove. We also summarize the mode frequencies at the Brillouin zone center (g = 0) and 
the zone edge {q = ir/a). As can be seen from the table, the transverse modes for all 
adsorbates considered have completely negligible dispersion throughout the Brillouin zone. 
We also find that u!ti{q) < <^r2(?) which is a consequence of a specific shape of the potential 
experienced by adsorbates in the groove. This means that the potential is "softer" with 
respect to adsorbate displacements in the 2;-direction (see Fig. 2). 

IV. SPECIFIC HEAT OF DENSE ONE-DIMENSIONAL ADSORBATE PHASES 

Knowing the dispersion relations for the adsorbate phonons, we can now proceed to 
calculate the specific heat pertaining to a dense phase of the adsorbates in a groove. The 
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heat capacity at constant volume, Cy, for a collection of independent oscillators (phonons) 
is given by^^ 



where is the Boltzmann constant, T is the temperature, and index i counts the harmonic 
oscillators of frequencies cui. In our case, the index i is replaced by the phonon wave vector, 
q and the phonon branch index s {s = L,T1,T2). Relation (11) is thus rewritten as 



{Lj./2'k) J J^/^ dq. Here, is the length of the adsorbate chain (or groove) and the number 
of adsorbates within a particular groove is N — L^/a. It is easy to check that in the limit 
of high temperatures {ksT » hijJs{q), s — L,T1,T2), the specific heat in Eq. (12) reduces 
to its classical value, Cy = SNks- 

In Fig. 3 we plot the specific heat for Xe,Kr and Ar adsorbates obtained from Eq. (12) 
by using information on the adsorbate phonons summarized in Eq. (6) and Table I. As can 
be seen, at high temperatures all the three curves tend to the classical limit. 

To separate the influence of transverse modes, we plot in Fig. 4 the three different 
contributions to the specific heat of a chain of Xe atoms, i.e. in the sum in Eq. (12) we 
consider only one mode (either L,T1, or T2). As can be inferred from this figure, the 
low temperature behavior of the specific heat (T < 4 K) is completely determined by the 
longitudinal mode. At temperatures higher than about 4 K (for Xe), all the modes in the 
sum in Eq. (12) must be considered. Due to the fact that the transverse modes have 
higher frequencies for Kr and Ar adsorbates, their contribution to the specific heat becomes 
significant at somewhat higher temperatures. The specific heat at low temperatures, where 
only the L mode contributes, is proportional to the temperature and the following relation 
holds: 




(11) 




(12) 



where we have transformed the sum over wave vectors into an integral as J2q 



Cv ~ 2m5NkB 



ksT 



(13) 



fiooiiq = tt/o) 
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The linear dependence of specific heat on temperature is a consequence of effectively ID 
behavior of the chain at low temperatures and can be easily seen from Fig. 3 at tempera- 
tures lower than a; 4K. Equation (13) also implies that the measurement of specific heat at 
low temperatures would yield an information on the maximum frequency of the longitudinal 
mode (Brillouin zone edge frequency). As this frequency crucially depends on the interadsor- 
bate potential [Eq. (5) and Eq. (6)], these measurements could yield an information on the 
effective interaction between the adsorbate atoms which may be modified by the presence 
of the substrate material, as discussed in Sec. III. 

V. SUMMARY 

We have studied the vibrational properties of hnear dense phases of adsorbates in the 
grooves of the carbon nanotube bundles. The three modes characterizing the adsorbate 
vibrations have been identified and their frequencies have been calculated for Xe, Kr and 
Ar adsorbates. We have also calculated the specific heat corresponding to the hnear phase 
and identified the temperature regimes in which the transverse modes begin to significantly 
contribute to the specific heat. The presented results should serve as a useful guidance in 
the measurements of specific heat of gases adsorbed in the nanotube materials. ^'^ 
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FIGURES 

FIG. 1. A sketch of a small bundle of carbon nanotubes with the one-dimensional adsorbate 

phase in one of its grooves. The choice of the coordinate system is denoted. The adsorbate atoms 
are denoted by small dark circles. 

FIG. 2. Full lines: External potential for Xe atom in the groove, (a) Potential along the 
z-direction. (b) Potential along the y-direction. Both curves pass through the absolute minimum 
of the potential located at {y = 8.5 A, z = 6.38 A). Dashed lines: Harmonic approximation to the 
potential. 

FIG. 3. Specific heat of Xe (thick full line), Kr (thick dashed line) and Ar (thick dotted 
line) adsorbate chains in a groove as a function of the temperature. The thin full line repre- 
sents the linear low-temperature behavior of specific heat for Ar as predicted by Eq. (13), with 
hujiiq = tt/o) = 4.93 meV (see Table I.). 

FIG. 4. Separate contributions of three modes of Xe chain to the specific heat. Dashed line: L 
mode. Dash-dotted line: T2 mode. Dash-dash-dotted line: Tl mode. Full line: total specific heat 
of the Xe chain. 
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TABLES 





Xe 


Kr 


Ar 


a [A] 


4.59 


4.03 


3.81 


fyy [N/m] 


7.13 


6.50 


5.58 


[N/m] 


1.01 


3.19 


2.58 


a [N/m] 


-0.003 


-0.003 


-0.0024 


P [N/m] 


1.10 


1.10 


0.87 


nujL{q = 0) [meV] 











hujL{q = vr/a) [meV] 


3.06 


3.83 


4.93 


fi(jjTi{q = 0) [meV] 


3.899 


4.660 


6.253 


huTi{q = 1^/0) [meV] 


3.895 


4.655 


6.247 


fiujT2{q = 0) [mcV] 


2.924 


3.264 


4.251 


h,UT2{q = [meV] 


2.919 


3.258 


4.244 



TABLE I. Force constants and phonon frequencies of L,T\ and T2 modes (Brillouin zone 
center and Brillouin zone edge) for Xe, Kr and Ar adsorbates in groove positions of the carbon 
nanotube bundle. 
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